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Abstract 

A baseline solution for CFD Point 1 (Mach 24) in the STS-107 accident investigation was modified to 
include effects of: (1) holes through the leading edge into a vented cavity; and (2) a scarfed, conical nozzle directed 
toward the centerline of the vehicle from the forward, inboard corner of the landing gear door. The simulations were 
generated relatively quickly and early in the investigation because simplifications were made to the leading edge 
cavity geometry and an existing utility to merge scarfed nozzle grid domains with structured baseline external 
domains was implemented. These simplifications in the breach simulations enabled: (1) a very quick grid generation 
procedure; and (2) high fidelity corroboration of jet physics with internal surface impingements ensuing from a 
breach through the leading edge, fully coupled to the external shock layer flow at flight conditions. These 
simulations provided early evidence that the flow through a two-inch diameter (or larger) breach enters the cavity 
with significant retention of external flow directionality. A normal jet directed into the cavity was not an appropriate 
model for these conditions at CFD Point 1 (Mach 24). The breach diameters were of the same order or larger than 
the local, external boundary-layer thickness. High impingement heating and pressures on the downstream lip of the 
breach were computed. It is likely that hole shape would evolve as a slot cut in the direction of the external 
streamlines. In the case of the six-inch diameter breach the boundary layer is fully ingested. The intent of externally 
directed jet simulations in the second scenario was to approximately model aerodynamic effects of a relatively large 
internal wing pressure, fueled by combusting aluminum, which deforms the corner of the landing gear door and 
directs a jet across the windside surface. These jet interactions, in and of themselves, were not sufficiently large to 
explain observed aerodynamic behavior. 


Introduction 

The Langley Aerothermodynamic Upwind 
Relaxation Algorithm (LAURA) 1 ' 2 was applied to 
coupled external - internal flow simulations for the 
Columbia investigation. For the case of leading edge 
damage, the simulations were generated to provide 
relatively high fidelity analyses of mass and energy 
ingestion rates through a breach as a function of hole 
size. They provided details of jet orientation and 
structure emanating from the breach into a vented 
cavity. These results were used to initialize 
engineering analyses of mass and energy flow 
through the internal structure of the wing and to 
better define the severity of the environment 
associated with the impinging jet. For the case of an 
externally directed jet through a deflected landing 
gear bay door, the simulations defined aerodynamic 
coefficients of the interaction to ascertain if such 


interaction could explain aerodynamic characteristics 
recorded in flight. 

Previous simulations of STS 1,2 and 5 with 
LAURA to resolve questions of a pitching moment 
anomaly 3 (STS 1) and to validate heating 
predictions 4 ' 5 were published in the AIAA Journal of 
Spacecraft and Rockets in 1994. Some of this 
material was also presented in an ITAR restricted 
document with more detail and dimensions. 6 There is 
no experimental or flight validation of the flow 
through the breach or the externally directed jet. The 
simulations were generated relatively quickly and 
early in the investigation on the NAS Cray S V 1 
because simplifications were made to the leading 
edge cavity geometry. These simplifications enabled 
(1) a very quick grid generation procedure; and (2) 
high fidelity corroboration of jet physics with internal 
surface impingements ensuing from a breach through 
the leading edge, fully coupled to the external shock 
layer flow at flight conditions. 


1 Approved for public release; distribution is unlimited. 



Software 

LAURA is a high fidelity analysis tool, 
specialized for hypersonic re-entry physics, utilizing 
state-of-art algorithms for computational fluid 
dynamic (CFD) simulations. Key elements of 
LAURA include Roe’s averaging 7 and Yee’s Total 
Variation Diminishing (TVD) 8 formulation of 
second-order, inviscid flux. Yee’s TVD formulation 
has been found to be exceptionally robust and 
Courant-number-independent using point-implicit 
relaxation for hypersonic flow simulations. The TVD 
algorithm uses a non-linear, minmod function as a 
flux limiter that maintains second-order accuracy 
away from extrema but can admit limit cycles in the 
convergence process, particularly in the vicinity of 
captured shocks. This occurrence usually manifests 
itself as a stalling of convergence at a very low error 
norm, essentially a benign ringing in the solution at a 
level that has no impact on aerothermodynamic 
quantities. Viscous flux is computed using central 
differences. 

Simulations 

All of the current simulations specify 
laminar flow of air in chemical nonequilibrium (N O, 
N 2 , CL, NO) and thermal equilibrium. Only CFD 
point 1 is simulated for internal cavity flow 
simulations with freestream conditions: V, , = 7350.6 
m/s, = 3.9005 10' 5 kg/m 3 , T M = 217.3 K, a = 
40.1681 deg., and time = 13.50.53.0. The external 
flow simulations are implemented on the shared, 
baseline grid developed for the External 
Environments Team. The baseline grid has been 
modified to allow coupling of the external flow with 
flow through a breach in the wing entering a vented 
cavity and with externally directed flow through a 
scarfed nozzle at the landing gear bay door. 

Holes in Panel 6 

The following material is the first 
application of LAURA to flow in a vented cavity. 
Consequently, geometric complexity was 
incrementally added to better understand ensuing 
flow physics and effects of boundary condition 
specification. 

Initially, quarter inch deep holes with 
varying diameters were inserted into the baseline, 
external flow solution at Panel 6. (See Figure 1.) The 
boundary condition at the bottom of the hole used a 
zero-order extrapolation from interior points for 
density and velocity and a specified reservoir 
pressure equal to 2 p„ : (approximately 0.1 psf). The 


flow through this boundary evolved to be 
substantially supersonic for hole-sizes greater than 
two inches. The vent backpressure is sufficiently low 
to maintain at least sonic conditions through the vent 
and is characteristic of leeside levels. 

The two-inch-hole simulation was then 
expanded to include an internal cavity extruded from 
the backside of panel 6. The cavity was roughly 1ft x 
lft x 1ft but its shape evolved from an ad-hoc grid 
generation process. A two-inch vent hole was 
included at the far end of the cavity. Flow through the 
breach is fully coupled to the external flow and the 
vent boundary condition was transferred to its new 
location. 



Figure 1: Streamlines over RCC panels on the 
leading edge and the location of the circular breach 
through the wing in panel 6. 

Procedure 

An in-house utility to insert structured grids 
defining Reaction Control System (RCS) jets through 
the outer mold lines (OML) of hypersonic vehicles 
was modified to define a quarter-inch deep, circular 
hole through the leading edge of shuttle orbiter wing. 
The modified surface grid blends smoothly with the 
baseline external grid as shown in Figure 2. A detail 
of the surface grid defining the hole is shown in 
Figure 3. 

Taking the inserted, rectangular domain 
surrounding the breach shown in Figure 2 and 
creating a quarter-inch offset orthogonal to the OML 
produced the near wall of the cavity. This surface 
was extruded into the interior in a direction defined 
by the axis of the hole. An initial cell size taken from 
the external flow domain was applied. A constant 
growth factor ( 1 .20) grows the grid off the wall for 
32 cells and then an inverse factor is applied for 
another 32 cells to resolve a boundary layer on the 
opposite wall. Sidewall boundary-layer resolution is 
achieved with a similar procedure in which additional 



points are added in directions defined by an 
extrapolated coordinate line. A view of this extruded 
cavity is shown in Figure 4. 



Figure 2: Modified surface grid in the vicinity of the 
hole. 



Figure 3: Structured grid within hole, approximately 
20 cells deep and 60 cells across. 

Immediate vent (No internal cavity) 


These simulations provided early evidence that the 
flow through a two-inch diameter (or larger) breach 
enters the cavity with significant retention of external 
flow directionality. A normal jet into the cavity was 
not an appropriate model for these conditions at CFD 
Point 1 (Mach 24). The breach diameters were of the 
same order or larger than the local, external 
boundary-layer thickness. High impingement heating 
and pressures on the downstream lip of the breach are 
computed. It is likely that hole shape would evolve as 
a slot cut in the direction of the external streamlines 
as shown in Figure 1 . In the case of the six-inch 
diameter breach the boundary layer is fully ingested. 
Summary of results follow. 
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Figure 4: Pressure in cutting plane across breach and 
view of surface grid defining the vented cavity. 


Two-inch hole (Completed 03-28-03) 


The contour plot of total enthalpy in Figure 
5 shows the boundary edge as the transition from 
orange to red color. The total enthalpy above the 
boundary layer edge, outside of viscous dissipation 
effects, is a constant. The boundary layer thickness is 
roughly % of the hole diameter. The pressure field 
and streamlines entering the breach (Figure 6) 
indicate significant retention of streamwise 
momentum after processing by a Prandtl-Meyer 
expansion emanating from the sonic line in the 
boundary layer. Surface heating rates exceed 300 
W/cm 2 (Figure 7) and surface temperatures approach 
3000 K (Figure 8) at the downstream lip of the 
breach where impingement is strongest. Ingested 
mass and energy rates are 0.000689 kg/s and 9.92 
kW, respectively. Boundary conditions (radiative 
equilibrium, finite catalytic wall) and examination of 
near surface grid quality (Re ce n of order 10 at 
impingement) suggest this heating rate is a lower 
bound. More data on grid convergence tests will be 
provided in the Panel 8 results to follow. 

Four-inch hole (Completed 03-31-03) 


The contour plot of total enthalpy in Figure 
9 is repeated for the four-inch diameter hole. The 
boundary layer edge impinges on the downstream lip. 
The pressure field and streamlines entering the 
breach (Figure 10) indicate even more retention of 
streamwise momentum as compared to the previous 
case for the smaller hole. Surface heating rates again 
exceed 300 W/cm 2 (Figure 11) and surface 
temperatures exceed 3000 K (Figure 12) at the 
downstream lip of the breach where impingement is 





strongest. Ingested mass and energy rates are 0.00255 
kg/s and 46 kW, respectively. 


Figure 8: Surface temperature contours in vicinity of 
the two-inch hole in Panel 6. 


Figure 5: Total enthalpy in cut plane above hole 
showing partial ingestion of the boundary layer. 
(Blue jags are artifact of picture format conversion.) 


Figure 9: Total enthalpy in cut plane above four inch 
hole showing significant ingestion of the boundary 
layer. (Blue jags are artifact of picture format 
conversion.) 


Figure 6: View of streamlines entering breach 
showing significant retention of external streamwise 
momentum. 


Figure 10: View of streamlines entering four-inch 
breach showing significant retention of external 
streamwise momentum. 
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Figure 7: Surface-heating contours in the vicinity of 
the two-inch hole in Panel 6. 
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Figure 1 1 : Surface-heating contours in the vicinity of 
the four-inch hole in Panel 6. 
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Figure 12: Surface temperature contours in 
vicinity of the four-inch hole in Panel 6. 


Figure 13: Total enthalpy in cut plane above six inch 
hole showing complete ingestion of the boundary 
layer. (Blue jags are artifact of picture format 
conversion.) 



Six-inch hole (Completed 04-01-03) 

The contour plot of total enthalpy in Figure 
13 is repeated for the six-inch diameter hole. The 
boundary layer edge and some inviscid flow at the 
maximum freestream total enthalpy above it are 
completely ingested. The lip is impacted by 
streamlines carrying the total freestream enthalpy and 
a relatively large total pressure. The pressure field 
and streamlines entering the breach (Figure 14) show 
the same trends for retention of streamwise 
momentum as compared to the previous case for the 
four-inch hole. Surface heating rates again exceed 
300 W/cm 2 (Figure 15) and surface temperatures 
exceed 3000 K (Figure 16) at the downstream lip of 
the breach where impingement is strongest. Of the 
three hole sizes, the present case should exhibit the 
largest heating because it is seeing undissipated flow 
from above the boundary layer edge. Ingested mass 
and energy rates are 0.00566 kg/s and 117 kW, 
respectively. 


Figure 14: View of streamlines entering six-inch 
breach showing significant retention of external 
streamwise momentum. 
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Figure 15: Surface-heating contours in the vicinity of 
the six-inch hole in Panel 6. 







Figure 16: Surface temperature contours in vicinity of 
the six-inch hole in Panel 6. 


Vented internal cavity 


These simulations allowed the jet to develop 
naturally within the cavity and to view initial 
impingement and deflections of flow to other 
surfaces. Jet orientation is a strong function of hole 
depth to diameter ratio. Summary of results follow. 

Two-inch hole to vented box (Completed 04-10-03) 

Figure 17 shows streamlines entering the 
breach and the ensuing reflected wave pattern 
producing compressions and expansions in the jet 
core. The jet is non-orthogonal to the back plane of 
the breach. It has a significant directional component 
defined by the external flow direction. An expansion 
off the upstream side of the hole overtakes a shock 
off the downstream side of the hole. The shock forms 
in response to supersonic flow from the upper portion 
of the boundary layer and the boundary layer edge 
impinging on the downstream lip. The Mach number 
distribution in Figure 18 provides additional detail of 
the jet entering the cavity. A sonic line sets up behind 
the inner edge of the breach lip shock indicating a 
choked condition. The interesting feature here is that 
the sonic line does not span the back plane of the hole 
but rather spans a smaller area defined by conver- 
gence of the streamlines entering from the circular 
border. Pressure coefficients in Figure 19 indicate the 
ambient pressure level in the cavity of approximately 
(1/8) p«= V„ 2 . This level is attained at the backplane 
of the hole. The dependence of cavity ambient 
pressure on vent hole size has not been explored. It is 
interesting to note that a specification of the vent 
pressure immediately behind the back plane of the 
hole results in approximately the same ambient 
pressure level in the cell preceding the outflow 
boundary. Pressures exceeding the external pressure 


of approximately (1/2) V TC 2 (C p = 1) would 
effectively block ingestion of external flow. 




Figure 17: Streamlines entering the breach and 
pressure levels in the cutting plane. External flow 
direction is from the right to the left. 



Figure 18: Mach number contours across the breach 
and in the jet entering the cavity. A sonic, choked 
condition sets up behind the lip shock over the 
downstream lip of the breach. 
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Figure 19: Pressure coefficient in vicinity of the 
breach. 





The temperature in the shock over the 
downstream lip approaches 8000 K in Figure 20. The 
corresponding total enthalpy profile is shown in 
Figure 21 that corresponds to the earlier result in 
Figure 5 except now the cavity beneath the hole is 
included. Conditions just below the boundary layer 
edge impact the lip to drive the large surface heating 
rates and surface temperatures (Figure 22). The 
impingement heating on the edge will tend to cut a 
slot along the surface in the external streamline 
direction. 

The oxygen content of ingested flow is 
almost fully dissociated as indicated in Figure 23. A 
strong catalytic heating effect would be expected on 
metallic cavity walls but the finite catalytic boundary 
condition on the external surfaces was continued into 
the cavity. Also, the radiative equilibrium wall 
boundary condition assumes cavity walls can radiate 
to free space whereas the actual closed environment 
would restrict radiative cooling. Conductive cooling 
through metallic walls is not considered. Given these 
conditions (expected to produce lower heating rates 
than actually encountered in flight) the surface 
temperatures on the cavity walls (Figure 24) already 
exceed the melting point of pure aluminum at the 
impingement zone. Note that the breach hole through 
the lower surface is not visible in this figure. The 
vent hole on the top surface still sees temperatures 
that would melt aluminum. 



Figure 20: Temperature in vicinity of the breach. 



Figure 21: Total enthalpy in vicinity of breach 
indicating level of ingestion of the external boundary 
layer. 
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Figure 22: Surface temperatures in vicinity of 
downstream lip of two-inch hole over cavity. 



Figure 23: Atomic oxygen mass fraction ingested 
through breach. 
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Figure 24: Surface temperature on cavity walls with 
identical cooling mechanisms as used for external 
flow. 


1/4-inch hole to vented box (Completed 04-28-03) 

A quarter-inch diameter hole was simulated 
to provide reference conditions in which the hole 
diameter was not larger than the hole depth and all 
hole dimensions were smaller than the local boundary 
layer thickness. In this case, the jet develops in an 
orthogonal direction to the wall as seen in the Mach 
number contour of Figure 25. Impingement heating 
and temperatures on the downstream lip are much 
more benign (Figure 26) because the energy content 
and total pressure of the impinging streamlines 
emerging from deep within the boundary layer is 
significantly dissipated. Temperatures entering the 
cavity (Figure 27) exceed 1000 K. 



Figure 25: Mach number contours of jet entering 
cavity through quarter-inch diameter hole. 



Figure 26: Surface temperatures in vicinity of 
quarter-inch hole. 



Figure 27: Jet temperatures in cavity bled from 
bottom of external boundary layer. 

Ten-inch Hole in Panel 8 


These simulations were initiated in mid-May 
to complement ongoing simulation work on higher 
fidelity internal grids and to insure that sufficient 
corroborating data for engineering analyses could be 
made available as quickly as possible. However, 
greater fidelity to the internal geometry was 
requested (as compared to the ad-hoc cavity) so that 
jet development lengths and splash to adjoining 
internal surfaces may be better simulated. 

A simulation plan was developed to 
maximize productivity in the limited, remaining time. 
First, the simulation of a ten-inch hole in Panel 8 into 
the ad-hoc cavity as described in the previous section 
was produced. (The grid in this case had to be post- 
processed by the VGM utility to compensate for 
problems associated with the larger ratio of local 
effective radius to hole diameter.) This simple 
simulation would next be used to initialize the cavity 
flow in a more realistic “chunnel” geometry 
generated concurrently. A grid sequencing process is 
used to complete the simulation and significant 





sharpening of the jet was observed in the transition 
from moderate to fine grid. When the jet boundaries 
sharpened, it was noted that the impingement 
location was under-resolved. An additional grid- 
processing step was required to bring a factor ten 
finer resolution on the wall at impingement. 

The chunnel geometry was built off of the 
ad-hoc cavity and expanded to include the domain 
behind panels 7-9. The back plane was defined by 
the actual spar location. Two vent slots with 
combined area of 66 in 2 were included on the leeside, 
back surface. An internal core system of blocks was 
approximately bounded by the inner extent of ribs at 
the panel interfaces to accommodate possible future 
extensions of internal complexity. A combination of 
an O- and C- grid topology was used to accommodate 
the complex internal structure. The C-grid provides 
the ability to capture the rib and spar geometry while 
offering a decoupling of the grid requirements for 
stretching and clustering at the internal wall surfaces. 
As shown in Figure 28 and Figure 29, the C-grid 
easily captures the necessary chunnel features. 

Within the C-grid, in the bottom block, the hole in 
the RCC is accommodated by simply replacing the 
block of the C-grid with an “embedded" O-grid, 
which enables accurate modeling of the flow entering 
the chunnel. The O- and C-grid junction enables the 
accurate modeling of the existing geometry while 
offering the expandability to increasing geometric 
complexity without restructuring and regenerating the 
entire chunnel volume grid. 

In the final simulation, the top vents were 
closed and the side core vents (combined area 79.3 
in 2 ) were opened to enable the splash from an 
impinging jet to spill over a rib and through a vent. 
Summary of results follow. 



Figure 28: Internal block structure merging ad hoc 
cavity and chunnel showing relative vent positions. 



Figure 29: Global view of chunnel block structure. 
Vented internal box (Completed 05-19-03) 

Mach number and streamlines entering the 
box through the ten-inch hole in panel 8 are shown in 
Figure 30. Streamlines exit a ten-inch diameter vent 
at the top of the box. The flow is strongly supersonic 
entering this cavity. The most interesting details of 
this simulation (that remain true for the subsequent 
chunnel simulations) are the high temperatures and 
pressures in the shock layer over the downstream lip 
of the hole. The high temperature (approaching 
12000 K in Figure 31) is easily explained because the 
lip is bathed in flow that carries the full freestream 
total enthalpy. The local pressure coefficient (C p = 

10) is higher than the stagnation point because the 
local streamlines here were processed by an oblique 
shock so that total pressure losses are less than those 
experienced in crossing a normal shock. Computed 
heating rates exceed 800 W/cm 2 in this case. When 
considering the factor of increase in pressure and 
factor decrease of local radius of curvature relative to 
the stagnation point one may expect this rate 
(ignoring ablative cooling) is a lower bound. 

Chunnel with slotted leeside vents (Completed 06- 
01-03) 

Results in this section discuss the solution 
after adaptation of the finest grid to concentrate more 
mesh to the jet impingement point. Error norms had 
dropped to order 0.4 (a five to six order of magnitude 
decrease from the initial solution off the coarse grid). 
Some unsteadiness is evident when monitoring the 
solution but there was insufficient time to gather 
meaningful statistics on frequency. The jet 
impingement point has moved very little over about 
10000 relaxation steps and it is believed that the basic 
flow physics presented here are correct. 



Figure 30: Mach number and streamlines in cavity 
behind panel 8 with ten-inch diameter hole and ten 
inch diameter vent. 


Figure 33: Mach number contours focused on 
supersonic domain. 



Figure 3 1 : Details of temperature in shock layer over 
downstream lip with ten-inch breach through panel 8. 
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Figure 34: Temperatures in cutting plane through 
chunnel showing high temperatures at impingement 
of jet on the side wall. 



Figure 32: Total enthalpy entering chunnel through 
Panel 8. Venting is through slots on top of chunnel. 


Figure 35: Comparison of chunnel surface pressures 
from medium grid and unadapted fine grid. 




Figure 36: Comparison of chunnel surface 
heating from medium grid and unadapted fine grid. 



Figure 37: Comparison of chunnel surface 
heating from unadapted fine grid and adapted fine 
grid on impingement surface. 

Details of the jet in the cutting plane through 
the chunnel are presented in Figure 32 (total 
enthalpy). Figure 33 (Mach number), and Figure 34 
(temperature). A tongue of high enthalpy flow 
persists up to the impingement point. A high 
temperature shock layer forms above this location. 
The jet initially scrapes the surface and then rises 
above it. 

A "strip" of high pressure (Figure 35) and 
high heating rate (Figure 36) is evident at the 
relatively narrow impingement location. These 
figures show the sharpening of the profile associated 
with refinement from the medium grid to the fine grid 
(a factor of 2 in the i, j, and k directions). Review of 
solution metrics at the impingement showed high cell 
Reynolds number and large jump between imposed 
wall temperature (500 K) and boundary cell center 
temperature (2880 K). Consequently, a grid 
adaptation was implemented to provide a target factor 
of 10 finer grid at the wall (i-direction) using a 
greater stretching factor for the same number of cells 
in the near wall blocks. Even with this refinement. 


the worst-case cell center temperature next to the 
wall was (1160 K). The convergence of peak 
impingement heating on these three sequential cases 
(not located at exactly the same surface mesh point) 
progressed from 116 W/cm 2 (medium grid) to 301 
W/cm 2 (fine grid) to 721 W/cm 2 (adapted fine grid. 
Figure 37). 

Chunnel with side vents (Completed 06-06-03) 

This side vent case (top vents now closed) 
was initialized from the previous top vent case with 
the finest grid and adaptation of grid on the 
impingement boundary. 

The error norm for this case jumped about 4 
orders of magnitude from the previous converged 
solution and returned to order 1 values at this 
"snapshot". An attempt to force the solution to 
advance more quickly with a constant, large time step 
eventually diverged. A smaller constant time step 
was again applied and the solution began to recover 
but insufficient time remained to drive it to 
convergence. The solution posted here is thought to 
have a jet impingement point that will continue to 
rise toward the side vent. The magnitude of jet 
oscillations cannot be estimated from the solution 
generated to date. A representative jet structure is 
presented in Figure 38 (total enthalpy) and Figure 39 
(pressure). The basic structure is very similar to that 
observed in the previous case. The expansion process 
to sonic conditions at the side vent is displayed in 
Figure 40 (pressure) and Figure 41 (Mach number) in 
the impingement zone. 

Top vent versus side vent internal surface 
pressures are compared in Figure 42. The equivalent 
heating comparison is presented in Figure 43 and 
Figure 44. Heating at impingement on the finest grid 
has remained between 500 and 700 W/cm 2 . 



Figure 38: Total enthalpy in chunnel with side vents. 




Figure 39: Pressures in chunnel with side vents. 



Figure 40: Detail of pressure at impingement location 
below side vent with supporting grid in the cut plane. 



Figure 41: Magnified view of Mach number 
expanding to sonic condition at side vent. 


Figure 42: Comparison of chunnel surface pressures 
with top and side vents. 



Figure 43: Comparison of chunnel surface heating 
with top and side vents. 



Figure 44: Logarithmic scaling of surface heating. 


Externally Directed Jet 

A baseline solution for CFD Point 1 (Mach 
24) in the STS-107 accident investigation was 
modified to include effects of a scarfed, conical 
nozzle directed toward the centerline of the vehicle 





from the forward, inboard comer of the landing gear 
door. The intent of the simulation is to approximately 
model effects of a relatively large internal wing 
pressure, fueled by combusting aluminum, which 
deforms the corner of the landing gear door and 
directs a jet across the windside surface. The present 
simulation does not include details of aluminum 
combustion but does consider extremes of internal 
conditions (p mtemal = 2 p„V«, 2 and T intemal = 4000 K) 
that are expected to provide an upper limit on this 
potential effect. (The associated analysis of flow 
through a breach in the leading edge of Panel 6 into a 
vented cavity indicates internal pressures of 
( 1/8) V. 2 fed by an external pressure of 

(1/2) V„ 2 ). The analysis includes interaction of the 

shock layer flow with the jet. The scarfed nozzle has 
a 3-inch (.0762m) diameter throat and a 1.3 ft 2 (.122 
m 2 ) elliptical footprint on the windside surface. 
Perturbations to baseline aerodynamic coefficients 
are expected to scale with throat area for the specified 
internal conditions because the interaction disturbs 
only a small region (roughly double the exit area) in 
the vicinity of the exit plane. Aerodynamic 
coefficient perturbations, provided in Table I, are of 
order 10" 3 to 10' 5 of right-half-body values. The 
exiting jet is substantially entrained in the windside 
boundary layer with relatively weak perturbations to 
the external flow. 


Table I Aerodynamic Coefficients 



Full- 

vehicle 

baseline 

A 

interactio 

n 

(new- 

baseline) 

A thrust 

A net 

axial 

.078175 

1.23 10 _/ 

3.07 10 _/ 

4.30 10' 7 

side 

0 

3.66 10° 

-8.12 lO' 3 

-4.46 10 -3 

norm. 

1.13226 

-3.03 lO' 3 

4.08 10 -5 

1.05 10' 3 

roll 

0 

-2.79 lO -3 

2.76 10 -5 

-0.03 10' 3 

pitch 

.027989 

3.98 lO' 3 

2.67 10 -6 

4.25 10 -5 

yaw 

0 

8.50 10' b 

-5.53 10' 6 

2.97 10' 6 


A force coefficient of 1.00 10" 5 corresponds to a force 
of 2.63 N (0.59 lb) for CFD Point 1 (Mach 24) 
conditions. A moment coefficient of 1.00 10" 5 
corresponds to a torque of 31.75 N-m (23.42 ft-lb). 

Jet Simulation Details 

Completed (04-23-03) 

A utility to carve out a structured subset of 
the baseline grid and replace it with grid that defines 
a scarfed nozzle has been applied. Several views of 
the modified grid are presented in Figure 45 to Figure 
47. Normally, additional grid refinement would be 


applied to provide smoother variation of cell size and 
better quality solution. In the present case, the 
original grid was employed to save time and because 
the grid was judged acceptable for the initial scoping 
of this problem. However, the severe change in 
resolution past the exit of the nozzle where the mesh 
width exceeds the nozzle width will completely 
dissipate solution profiles once the jet crosses to the 
coarse mesh. Because the jet never significantly 
pierced the boundary layer even where the flow was 
highly resolved (see streamlines from core of jet in 
Figure 48) a subjective judgment was made that 
additional refinement was unlikely to produce 
significant additional changes and resources should 
be focused on other scenarios. Note the analysis was 
conducted on the right side of the vehicle but all 
aerodynamic perturbations in Table I have been 
transformed for a jet on the left wing. The center of 
the jet exit plane on the windside surface is located at 
(x,y,z) = (26.54m, 2.91m, 6.90m). This location 
places the far end of the exit near the forward, 
inboard corner of the landing gear door. The 
simulation does not include the effects of a deformed 
door extending into the boundary layer. However, the 
jet itself does displace the external flow in a similar 
manner. 



Figure 45: Global view of embedded nozzle grid 
within the baseline grid. 




Figure 46: View of surface mesh in vicinity of 
scarfed nozzle exit. 



Figure 47: View of grid defining the inner surface of 
the nozzle and the feeder block. 


The inflow conditions are computed 
assuming isentropic expansion from a user specified 
stagnation pressure (2p_,_ V..,. 2 ) and stagnation 
temperature (4000 K). The gas is specified as 
dissociated air with molecular nitrogen mass fraction 
equal to 0.77 and atomic oxygen mass fraction equal 
to 0.23. (Conditions in a vented cavity fed by a 
breach through RCC panel 6 are later used to provide 
context for these assumptions.) Inflow conditions at 
the throat are held constant and feed the nozzle 
blocks. Boundary values at the center of the jet at the 
interface between the feeder blocks and nozzle blocks 
(See Figure 48) are defined in Table II. 


Table II - Jet Inflow Conditions Near Throat 



dimensional 

non- 

dimensional 

reference 

p 

1.591 10" 3 kg/m 3 

40.79 

= 3.9 1 0' 5 kg/m 3 

u 

-0.484 m/s 

-6.58 10' 5 

V 00 = 7350.6 m/s 

V 

1325.82 m/s 

0.1804 

V M 

w 

-525.323 m/s 

-0.0715 

V„ 

p 

1938.15N/m 2 

0.9198 

P°° Vco 2 

T 

3497 K 



H 

8.95 MJ/kg 

.1656 

vj 

M 


1.11 



Z 



Figure 48: Streamlines in jet core emanating from 
feeder block at nozzle throat. Streamlines are 
substantially entrained in the boundary layer. 
Pressures contours are in cutting plane through 
middle of nozzle. See Figure 49 for global 
orientation. 

The extent of any significant interaction is 
confined to the near exit region; it was not necessary 
to include the opposite side of the vehicle in the 
simulation for the given reservoir (landing gear bay) 
conditions. Streamlines from the jet are turned 
downstream well outboard of the plane of symmetry 
as seen in Figure 49. The streamlines would likely 
penetrate further toward the centerline if grid quality 
were improved. A modest push of the plane of 
streamlines within the open part of the nozzle in the 
external flow direction is evident in Figure 50. As 
each streamline rises above the level of the outer 
mold line it is quickly entrained in the external 
boundary layer. Higher pressures upstream of the 
nozzle are associated with the compressive action of 
the boundary layer being lifted by the jet. Lower 
pressures persist somewhat further downstream in the 
wake of the jet. The negative normal force 
aerodynamic coefficient from only interaction effects 






(pressures and shear forces on solid surfaces) in 
Table I indicates the integrated wake effect is slightly 
higher than the upstream compression effect. 



Figure 49: Global view of pressure in shock layer in 
cutting plane through nozzle and streamlines from 
nozzle flowing downstream. 


Figure 50: Overhead view of streamlines emanating 
from vertical line across jet at the feeder block. 
Pressure levels on the surface are colored from red 
(high) to blue (low). 

Summary 

A baseline solution for CFD Point 1 (Mach 
24) in the STS-107 accident investigation was 
modified to include effects of: (1) holes through the 
leading edge into a vented cavity; and (2) a scarfed, 
conical nozzle directed toward the centerline of the 
vehicle from the forward, inboard corner of the 


landing gear door. The simulations were generated 
relatively quickly and early in the investigation 
because simplifications were made to the leading 
edge cavity geometry and an existing utility to merge 
scarfed nozzle grid domains with structured baseline 
external domains was implemented. These 
simplifications in the breach simulations enabled: (1) 
a very quick grid generation procedure; and (2) high 
fidelity corroboration of jet physics with internal 
surface impingements ensuing from a breach through 
the leading edge, fully coupled to the external shock 
layer flow at flight conditions. 

These simulations provided early evidence 
that the flow through a two-inch diameter (or larger) 
breach enters the cavity with significant retention of 
external flow directionality. A normal jet directed 
into the cavity was not an appropriate model for these 
conditions at CFD Point 1 (Mach 24). The breach 
diameters were of the same order or larger than the 
local, external boundary-layer thickness. High 
impingement heating and pressures on the 
downstream lip of the breach were computed. It is 
likely that hole shape would evolve as a slot cut in 
the direction of the external streamlines. In the case 
of the six-inch diameter breach the boundary layer is 
fully ingested. The boundary-layer edge impinges on 
the downstream lip of a four-inch diameter hole. The 
boundary layer is not fully ingested for a two-inch 
diameter hole but significant retention of the external 
flow momentum is evident in the jet. Only the bottom 
part of the boundary layer is ingested through a 
quarter-inch diameter hole and the ensuing jet into 
the cavity is nearly orthogonal to the mold line. 

The intent of externally directed jet 
simulations was to approximately model 
aerodynamic effects of a relatively large internal 
wing pressure, fueled by combusting aluminum, 
which deforms the corner of the landing gear door 
and directs a jet across the windside surface. These 
jet interactions, in and of themselves, were not 
sufficiently large to explain observed aerodynamic 
behavior. The simulated geometry is thought to be a 
reasonable approximation to the speculated scenario. 
It seems unlikely that a significantly more effective 
nozzle shape would have evolved. 

The present simulation assumed a steady 
flow with chamber stagnation conditions sustained by 
combustion of aluminum (though the details of 
aluminum combustion chemistry are ignored). The 
results of the cavity simulation show that the 
assumed chamber pressures would choke off flow of 
oxygen from a breach to sustain the combustion. 
However, a pulse jet scenario is likely sustainable, 
even with higher stagnation pressures and 
temperatures in the chamber, for short, periodic 
bursts. The natural periodic frequency of such a 





process cannot be estimated here because of the 
complexity (and time dependency) of internal flow 
paths. 
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